Non-unitary triplet pairing in the centrosymmetric superconductor LaNiGa 2 



(N 

o 

£ 

(N 



o 
o 

I 

u 

Oh- 
o3 



c 

o 
o 



> 
in 
o 
<^ 
in 

O 

(N 



X 



A.D. Hillier, 1 J. Quintanilla, 1 ' 2 B. Mazidian, 1 * 3 J. F. Annett, 3 and R. Cywinski 4 

1 ISIS facility, STFC Rutherford Appleton Laboratory, 
Harwell Science and Innovation Campus, Oxfordshire, 0X11 OQX, UK 
2 SEPnet and Hubbard Theory Consortium, School of Physical Sciences, University of Kent, Canterbury CT2 7NH, UK 
H. H. Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol BS8 1TL, UK 
^School of Applied Sciences, University of Huddersfield, Queensgate, Huddersfield, HD1 3DH, UK 

(Dated: March 8, 2013) 

Muon spin rotation and relaxation experiments on the centrosymmetric intermetallic supercon- 
ductor LaNiGa2 are reported. The appearance of spontaneous magnetic fields coincides with the 
onset of superconductivity, implying that the superconducting state breaks time reversal symme- 
try, similarly to non-centrosymmetric LaNiC2. Only four triplet states are compatible with this 
observation, all of which are non-unitary triplets. This suggests that LaNiGa2 is the centrosym- 
metric analogue of LaNiC2. We argue that these materials are representatives of a new family of 
paramagnetic non-unitary superconductors. 



Symmetry breaking is a central concept of physics for 
which superconductivity provides one of the best un- 
derstood paradigms. In a conventional superconductor 
[l| gauge symmetry is broken, while unconventional su- 
perfluids and superconductors break other symmetries 
as well [3]. Examples include 3 He Q, cuprate high- 
temperature superconductors jij, the ruthenate Sr 2 Ru04 
[j| and more recently, non-centrosymmetric LaNiC2 0. 
The latter has weak spin-orbit coupling (SOC) Q, low 
symmetry and is a non-unitary superconductor. In a non- 
unitary superconductor the pairing states of the spin- 
up and spin-down Fermi surfaces are different. At the 
instability, a spin-up superfluid can coexist with spin- 
down Fermi liquid. While non-unitary triplet super- 
conductivity is well-established in ferromagnetic super- 
conductors its occurrence in parmagnetic LaNiC2 
remains puzzling. Here we provide experimental evi- 
dence of this phenomenon in another, compositionally 
related, but centrosymmetric superconductor: LaNiGa2. 
We also advance an explanation in terms of a coupling 
between triplet instabilities and paramagnetism that is 
quite generic and for which these two could provide the 
first examples of what might be a larger class of materi- 
als. 

In general unconventional pairing can be difficult to 
establish in any given material. However evidence for 
time-reversal symmetry (TRS) breaking in particular can 
be shown through the detection of spontaneous but very 
small internal fields Q. Muon spin relaxation/rotation 
(/iSR) is especially sensitive for detecting small changes 
in internal fields and can easily measure fields of 0.1 G 
which corresponds to wO.01 hb- This makes /iSR an 
extremely powerful technique for measuring the effects 
of TRS breaking in exotic superconductors. Direct ob- 
servation of TRS breaking states is extremely rare and 
spontaneous fields have been observed in this way only in 
a few systems: PrOs 4 Sbi 2 @, Sr 2 Ru0 4 (where TRS 
breaking was subsequently confirmed by optical measure- 
ments [HI), B-phase of UPt3 [12] (although not without 



controversy 0,01), (U.Th)Bei 3 [H] and more recently 
LaNiC 2 @, PrPt 4 Gei2 03] and (PrLa)(OsRu) 4 Sbi 2 [HJ. 
For examples of other systems where the effect is not ob- 
served see Refs. (lS - 21 1 . Broken TRS in superconduc- 
tors is especially interesting, because it implies not just 
unconventional pairing, but the existence of two-fold or 
higher de gen eracy of the superconducting order parame- 
ter space [22|. 

The observation of broken TRS in LaNiC 2 was partic- 
ularly surprising because of the low symmetry of this or- 
thorhombic, non-centrosymmetric, materialp]. Symme- 
try analysis has shown that the low dimensionality of this 
structure, with C 2t , point group, gives rise to only 12 pos- 
sible gap functions. Of these only 4 break TRS and these 
are all non- unitary triplet pairing states [f|. These four 
gap functions are all derived from one dimensional irre- 
ducible representations of the point group, implying that 
the only possible order parameter degeneracy is derived 
from the triplet Cooper pair spin orientational degree of 
freedom. A subsequent analysis of the effects of SOC on 
this system 0] shows that the SOC always lifts this fi- 
nal degeneracy, leading to a completely non-degenerate 
order parameter space, which would not be expected to 
allow spontaneous breaking of TRS at the superconduct- 
ing transition temperature, T c [22]. The only way to 
reconcile this with the experimental observations of bro- 
ken TRS in this material is to assume that the effect of 
SOC is weak on the relevant electron states at the Fermi 
level in this material. Additional experimental evidence 
for unconventional pairing in LaNiC 2 has been reported 
recently [23| . 

In this letter we report ^tSR results on the centrosym- 
metric superconductor LaNiGa 2 showing that TRS is 
broken on entering the superconducting state. This 
is a centrosymmetric material, which crystallises in 
the NdNiGa 2 orthorhombic structure, with space group 
Cmmm (D 2 h) [3] (see Fig. Q]). Magnetisation and 
heat capacity measurements have previously shown that 
LaNiGa 2 is a paramagnetic superconductor, with a T c 
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FIG. 1: (color online) The orthorhombic crystal structure 
of LaNiGa2. The red spheres (largest) are La, blue spheres 
(smallest) are Ni and the black spheres (medium) are Ga. 



onset of 2.1 K[2J|. Heat capacity measurements have 
shown a specific heat jump AC/"/T c «1.3l[34j, which 
is slightly lower than the expected BCS value of 1.43, 
and the temperature dependence has the conventional 
BCS exponential form. Specific heat jumps in TRS- 
breaking superconductors are sometimes lower and some- 
times higher than the BCS value, for example Sr2Ru04 
LaNiC 2 (li| and PrPt 4 Gei 2 [27|], respectively. Below 
we analyse the possible order parameter symmetries. Al- 
though LaNiGa2 has a different point group to LaNiC 2 
and is centrosymmetric, it is similar to LaNiC 2 in that it 
has only 12 possible superconducting symmetries For 
both materials these are all derived from one-dimensional 
irreducible representations. In particular, only four of 
these states break TRS and these are all non-untiary 
triplet pairing states (see Table Hk). 

The sample was prepared by melting together stoichio- 
metric amounts of the constituent elements in a water- 
cooled argon arc furnace. The fiSK experiments were 
carried out using the MuSR spectrometer in longitudinal 
and transverse geometries. At the ISIS facility, a pulse 
of muons is produced every 20 ms and has a FWHM of 
^70 ns. These muons are implanted into the sample and 
decay with a half-life of 2.2 /is into a positron which is 
emitted preferentially in the direction of the muon spin 
axis and two neutrinos. These positrons are detected and 
time stamped in the detectors which are positioned either 
before, F, or after, B, the sample for longitudinal (relax- 
ation) experiments. The forward and backward detectors 
are each segmented into 32 detectors. Using these counts 
the asymmetry in the positron emission can be deter- 
mined and, therefore, the muon polarisation is measured 
as a function of time. For the transverse field experi- 
ments, the magnetic field was applied perpendicular to 
the initial muon spin direction and momentum. For a 
more detailed description of the different instrumental 

ilioj]. 



geometry can be found in Ref. 

The sample was powder mounted onto a 99.995+ % 
pure silver plate. Any muons stopped in silver give 
a time independent background for longitudinal (relax- 
ation) experiments. The sample holder and sample were 
mounted onto a TBT dilution refrigerator with a temper- 



ature range of 0.045-4 K. The stray fields at the sample 
position are cancelled to within 1 /iT by a flux-gate mag- 
netometer and an active compensation system controlling 
the three pairs of correction coils. The transverse field 
/iSR (TF-/iSR) experiment was conducted with applied 
fields between 5 mT and 60 mT, which ensured the sam- 
ple was in the mixed state. Each field was either applied 
above the superconducting transition and the sample was 
then cooled to base temperature (FC) or the sample was 
first cooled to base temperature and then the field was 
applied (ZFC). The sample was cooled to base temper- 
ature in zero field and the /iSR spectra were collected 
upon warming the sample while still in zero field. 

The MuSR spectrometer comprises 64 detectors. In 
software, the data is mapped to two orthogonal virtual 
detectors each characterised by a phase offset ip. The 
resulting 2 spectra were simultaneously fitted with a si- 
nusoidal oscillating function with Gaussian relaxation: 
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where the ith index denotes the sample and background 
contributions, respectively, Ai is the initial asymmetry, 
<Ji is the Gaussian relaxation rate and is the muon 
spin precessional frequency. The background term comes 
from those muons which were implanted into the silver 
sample holder and therefore this oscillating term has no 
depolarisation, i.e. <r 2 =0.0 ^s _1 , as silver has a negligible 
nuclear moment. Fig. [2] shows a typical spectrum for 
LaNiGa 2 with an applied field of 40 mT at 50 mK after 
being FC. 

As the muon spin rotation arising from the field distri- 
butions associated with the flux line lattice is indepen- 
dent of that arising from the nuclear moments we can 
write af = a 2 ac + er 2 . o n is assumed to be constant in 
this temperature region, and is determined from mea- 
surements just above T c . Each data point was collected 
after field cooling the sample from above T c . The field 
dependence of <j sc (in /is) is related to the superconduct- 
ing penetration depth A (in nm) and coherence length £ 
via the relation 



a sc = 4.83 x 10 4 (1 - 6)(1 + 1.21(1 - Vbf)X' 



(2) 



where b = B/Bc2 is the ratio of applied field to upper 
critical field. From this we have determined A and Bc 2 
and hence £ to be 350(10) nm, 410(3) mT and 28(3) nm 
respectively (see Fig. [3]). This shows that LaNiGa 2 is a 
type II superconductor, with a superconducting electron 
density and effective superelectron mass of 9 x 10 26 m~ 3 
and 3.9 m e , respectively. More details on these calcula- 
tions can be found in Ref [3l|. Now let us consider the 
longitudinal /j,SK data. The absence of a precessional 
signal in the /iSR spectra at all temperatures confirms 
that there are no spontaneous coherent internal mag- 
netic fields associated with long range magnetic order in 
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FIG. 2: (Color online) The upper graph is a typical muon 
asymmetry spectra in LaNiGa2 taken in a transverse field of 
40 mT at 0.05 K (shown in the rotating reference frame (RRF) 
of 6.0 MHz. The line is a fit to the data using Eqn. [TJ For 
clarity, only one of the two virtual detectors have been shown. 
The lower graph is the zero field /iSR spectra for LaNiGa2. 
The blue symbols are the data collected at 56 mK and the 
red symbols are the data collected at 3.0 K. The lines are a 
least squares fit to the data. 



LaNiGa2 at any temperature. In the absence of atomic 
moments muon spin relaxation is expected to arise en- 
tirely from the local fields associated with the nuclear 
moments. These nuclear spins are static, on the time 
scale of the muon precession, and are randomly orien- 
tated. The depolarisation function, G z ( t), can be de- 
scribed by the Kubo-Toyabe function [32| 
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where (7/7^ is the local field distribution width and 7^ = 
2ir x 135.5 MHz T _1 is the muon gyromagnetic ratio. The 
spectra that we observed for LaNiGa2 are well described 
by the function 



G z (t) = A G* T {t) exp(-At) + A bckgrd , 



(4) 



where Aq is the initial asymmetry, Abckgrd is the back- 
ground, and A is the electronic relaxation rate (see Fig. 
[2]). It is assumed that the exponential factor involving 
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FIG. 3: (Color online) The field dependence of a at 50mK, 
after being field cooled. The line is a fit to the data using 
Eqn. [2] 
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FIG. 4: (color online) The left graph shows the tempera- 
ture dependence of a, for LaNiGa2 in zero-field, which clearly 
shows the spontaneous fields appearing at T c =2.1 K (shown 
has the vertical line). The line is fit to the data using an 
approximation [3^1 to the BCS order parameter for a e . The 
right graph shows the temperature dependence of the elec- 
tronic relaxation rate, A, for LaNiGa2 in zero-field, which 
shows no temperature dependence. 



A arises from electronic moments which afford an en- 
tirely independent muon spin relaxation channel in real 
time. The only parameter that shows any temperature 
dependence is <r, which increases rapidly with decreasing 
temperature below T c (see Fig. H|. We interpret this in- 
crease in a as a signature of a coherent internal field with 
a very low frequency as discussed by Aoki et al. Q for 
PrOs4Sbi2- This increase in a has been modelled, assum- 
ing that there are uncorrelated, by t(T) 2 = a 2 n + a e (T) 2 , 
where er„ and <j e are the nuclear and electronic contri- 
butions respectively. The temperature dependence of er e 
agrees with the BCS order parameter (see Fig. @|. 

Tet us now discuss the implications of this result for 
the pairing symmetry. The group theoretical analysis 
for the D2/1 point group of this system has already been 
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in our system is of the form 



TABLE I: The upper table (a) shows the gap function of 
the homogeneous superconducting states allowed by symme- 
try, for weak spin-orbit coupling. We have used the stan- 
dard notation A(k) = A(k)ia y for singlet states and A(k) = 
i [d(k).<x] cj y for triplets, where <r = (a x , a y ,a z ) is the vector 
of Pauli matrices. Here X, Y and Z are any basis functions 
in the Brillouin zone of odd parity, such as sin k x , sin k y and 
sinfc z , respectively. Only the four non-unitary triplet states 
are compatible with our observation of broken TRS. The lower 
table (b) shows the homogeneous spin-triplet superconduct- 
ing states allowed by symmetry at T c , for the case of strong 
spin-orbit coupling. The spin singlet gap functions are the 
same as in the weak spin-orbit case. The notation of the pos- 
sible gap functions is the same as the upper table, except that 
here A, B and C are (real) constants determined by the mi- 
croscopic gap equation. Clearly none of these states breaks 
TRS at T c . 



investigated [J]. For the simplest case, where transla- 
tional symmetry is not broken and SOC does not play a 
role, this point group has a total of 8 irreducible repre- 
sentations. This leads to 12 possible order parameters, 
as given in Table [T^l. Of these 12, 8 are unitary and 4 
are non-unitary. Only the 4 non-unitary order parame- 
ters have a non-trivially complex order parameter that 
can break TRS. In the case where SOC is large there are 
only 4 possible states (see Table |TJd) and none of them 
break TRS. Therefore, like LaNiC2, LaNiGa2 must be a 
non-unitary triplet superconductor with weak SOC. As 
we have predicted for LaNiC2@ if the SOC is not zero 
then a split transition would be expected. 

Until the discovery of non-unitary triplet pairing in 
LaNiC2 this state had only been confirmed in ferromag- 
netic superconductors Q. The additional observation 
of non-unitary triplet pairing in LaNiGa2 brings up the 
question of how a triplet superconductor whose normal 
state is paramagnetic could favour this state. The usual 
Landau free energy describing a triplet pairing instability 



F = a\r 1 \ 2 + -\r 1 \ 4 + b'\r 1 xr 1 * 



where ry is the order parameter, which relates to the d 
vector through d (k) = rjT (k) [the possible functional 
forms of r (k) are given in Table H^,] . The triplet insta- 
bility takes place when a = 0, which determines T c and is 
independent of whether pairing is unitary or non-unitary. 
Below T c , the second of the quartic terms decides which 
of the two states is most stable. The criterion for non- 
unitary triplet pairing is Q 



b' < 0. 



(6) 



On the other hand for a paramagnet there must be an 
additional term coupling 77 to the magnetization m. On 
symmetry grounds the simplest form of the free energy 
that takes this into account is 



F 



a\v\ 



\r]f + b' \t) X T7*| 2 + b"m-(ir) x rj*) . 

(7) 

Here \ is the normal state susceptibility. 

For given r\ the last term on the right hand side 
of (JTJl describes an effective magnetic field h e // = 
—b" (ir] XT]*) coupled to m. This field vanishes for uni- 
tary triplet pairing, but in the non-unitary case it induces 
a magnetization 



-Xb" (irj xr]*). 

a/2 



(8) 



Below T c , r] ~ (T c - T) i/Z whence m - T c -T. This 
subdominant order parameter lowers the energy of the 
non-unitary state compared to the unitary one. Indeed 
substituting ([H} into (J7|) we recover the simpler expression 
([5|) but with the 6' coefficient replaced with b' — 6" 2 x/2. 
The condition ^ then becomes 



b" 2 X /2 < 0. 



(9) 



For a paramagnet the second term on the left hand side 
is always negative, favouring non-unitary triplet pairing 
states. This effect would be expected to be strongest 
in proximity to a Stoner instability. We note that in 
superconducting ferromagnets @ the same coupling term 
exists and stabilizes non-unitary triplet pairing states by 
increasing their T c relative to unitary states. 

In conclusion, zero field and transverse field /iSR exper- 
iments have been carried out on LaNiGa2 . The zero field 
measurements show a spontaneous field appearing at the 
superconducting transition temperature. This provides 
convincing evidence that time reversal symmetry is bro- 
ken in the superconducting state of LaNiGa2. Symme- 
try analysis implies non-unitary triplet pairing, in close 
analogy with the non-centrosymmetric superconductor 
LaNiC2- We propose that these materials could repre- 
sent a new class of superconductors where a triplet su- 
perconducting instability of a paramagnetic state gives 
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rise to non-unitary pairing through a generic coupling to 
the magnetization. 
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